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1. Introduction 

The photoproduction reaction of neutral pions on the proton is a reaction which shows a strong 
effect due to isospin breaking. The electric multipole Eo+ exhibits an exceptionally strong cusp at 
the 7i~^n threshold (see for instance Ref. [1]). The strength of this cusp is intimately related to the 
charge exchange scattering length of pion-nucleon scattering. Therefore, accurate experimental 
data of the photoproduction reaction allows one to access the pion-nucleon scattering lengths. 

Along the lines of previous work which analyzed the cusps in A' ^ 371 and rj' r\nn de- 
cays [2-4], we construct a nonrelativistic theory which provides a rigorous framework to describe 
the structure of the cusp order by order in a perturbative expansion in terms of small momenta 
and effective range parameters of pion-nucleon scattering as well as threshold parameters of the 
photoproduction reaction [5]. By construction, the theory correctly reproduces the low-energy 
singularities in the Mandelstam plane. 

The cusp in neutral pion photoproduction has been studied before. Ref. [6] introduces a two- 
parameter model, which captures the most important leading effect of the cusp. In Ref. [7], a 
coupled channel 5-matrix approach is used to investigate the cusp structure. 

2. Multipole decomposition 

Some basic relations and definitions used in the analysis of pion photoproduction are collected. 
We calculate the matrix element for the process p{p\) + y{k) p{p2) + T^^iq) at leading order in 
the electromagnetic coupling e, 

{p2,qoxA\py,km) = -i{2nf 8^^\Pf - Pi)u{p2,t')ei,J^u{put) , (2.1) 

where Pi and Pf denote the total four momentum in the initial and in the final state, respectively 
and stands for the polarization vector of the photon. 

To analyze the photoproduction reaction of pions, electric and magnetic multipoles are usually 
introduced. To this end, the ampUtude is written in terms of two component spinors t,t and Pauli 
matrices [8], 

5" = /T-e?"! + T-qT - (k X e)3'2 + 'T-kq-e J3 + /T-qq-£?'4. (2.2) 

The hat denotes unit vectors. The 9", are decomposed into electric and magnetic multipoles with 
the help of derivatives of the Legendre polynomials Pi{z) [8], 

3^1 = £ [/M,+ + Ei^P'i^, (z) + [(/ + mi- + Ei-]P'i-M , 
1=0 

3'2 = l,[{l + mi+ + lMi-]P'i{z), 
1=1 

3^3 = [£/+ -M/+]P/;i(z) + [Ei_ +Mi_]P'/_i{z) , 
1=1 

% = Y. [Mi+ - Ei+ - Ml. - EiJ\P'i'{z) . (2.3) 
/=i 

The discussion is restrained to the center of mass frame in the rest of the article. 
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3. Nonrelativistic framework 

To describe the behavior of the multipoles close to threshold - where the energy of the pro- 
duced pion and of the proton are small - a nonrelativistic calculation is justified. Furthermore, it 
offers the advantage that all the masses can be set to their physical value. Therefore, all the poles 
and branch points appear at the correct place in the Mandelstam plane. Moreover, the interaction 
of the nucleon and the pion is described by effective range parameters, which allows one to directly 
access the pion-nucleon scattering lengths. 

The covariant formulation of nonrelativistic field theories introduced in 
Refs. [2-4] is used here since it incorporates the correct relativistic dispersion law for the parti- 
cles. The nonrelativistic proton, neutron and pion fields are denoted by Xj/, % and iZk, respectively. 
The kinetic part of the Lagrangian after minimal substitution takes the form (see Ref. [9]) 

^kin = £(/4AW±;r± - i{D,^N±n±)^n± - 24w|;r±) 

+ 2x^Wn{idt-Wn)x + 2%lWQ{id,-Wo)%Q , (3.1) 



Wo = ^M^o-A, Wn = yJml-A, D,n± = {dt t ieAo)n± , 



with 



Dn± = (V ± ieA)n± , Dy/- = (V + ieA)^- (3.2) 

Note that since the photon is treated as an external field, its kinetic term is absent. 

4. Power counting 

Close to threshold, the momenta of the incoming proton and photon are of the order of the 
pion mass whereas the outgoing particles have very small momenta. Therefore, we count momenta 

of the outgoing pion and the outgoing proton as a small quantity of 0(e) and the momenta of 
the incoming proton and of the photon as 0(1). All the masses are counted as 0(1). The mass 
differences of the charged and neutral pion, A„ = — and of the proton and the neutron, 
Af^ = m^ — m^p are counted as 0{e^). At first sight, this counting scheme seems to lead to infinitely 
many terms already in the leading order p + Y^ p + Ti'^ Lagrangian =Sfy because derivatives on 
the incoming fields are not suppressed. However, since the the modulus of the momentum of the 
incoming particles, |k|, can be expanded in the small momentum |q|, 

to = ^^^, t, = |^^, y = ^, (4.1) 

„ 2 1-l-y 4M^{l+y) nip 

one obtains a valid power counting scheme. The Feynman rule in momentum space of every 
operator of order e*^ with a given arbitrary number of derivatives can be expanded in powers of 
the small momentum q, yielding one term of order e*^ without any momenta of the incoming fields 
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present and subsequent higher order terms. Doing this for every operator of 0{e^), all the resulting 
leading order terms without any momentum dependence can be described by one operator of order 
0{e'^) in the interaction Lagrangian. The same procedure leads to finite numbers of operators at 
any given higher order in e. The derivatives on the incoming fields are only needed to generate unit 
vectors in the direction of the incoming photon. This shows that the nonrelativistic theory is not 
capable of predicting the dependence on |k| even at threshold. 

An additional generic parameter a is introduced to count the pion-nucleon scattering vertices. Every 
pion-nucleon interaction vertex counts as a quantity of order 0{a) since the coupling constants are 
proportional to the pion-nucleon scattering threshold parameters, which are small. The perturbative 
expansion is therefore a combined expansion in e and a. 

5. Interaction Lagrangian 

The Lagrangian needed for the calculation of the amplitudes for pion photoproduction reads 
.if = ^kin + + ■^TtN^ where .if^,„ denotes the kinetic part, incorporates the interaction with 
the photon field and and S^TtN describes the pion-nucleon sector. 

In the pion nucleon sector, the leading terms of the Lagrangian have been given before in 
Ref. [10]. First, some notation is introduced in order to write the Lagrangian in a compact form. 
For every channel n, we collect the charges of the outgoing and the incoming pions in the vari- 
ables v and w, («;v,w): (0;0,0),(1;0,+), (2;+,+), (3;0,0), (4;-,0), (5; -,-), thereby assigning 
unique values to the variables v and w once n is given. The Lagrangian reads 



7(0.5} 7(1,4} 




V'(l,4} ^{2,3} 



(5.1) 



with the abbreviation f Ag = fAg + {Af)g. 

For the photon is treated as an external vector field A which is odd under parity and time- 
reversal transformations. One obtains for the gauge invariant Lagrangian 

££!f^ = -iGf V^'!^VV''''£*^V^''4 - iGf\if^%^\ifE^ATtl 

+ /Gi|V"^T'"TVB'AV'"7r,j + iGf^Y^^'"^'¥'^'"B'v'"'"7il 

- /gS3V^tVv-''£*^AV*^4 - /Gi4V^TVv^'''"£^V'=''";ro^ 

- iGfj \\fUj\\f¥E''Vj''^nl . (5.2) 
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The upper index on the couphng constants is introduced for later convenience. Here, the notation 
yim-ik = y'iy!2 . . . y;* jg used. Since the structure of the Lagrangian for the other required channel 
PY — > nn^ stays the same, one only has to replace the coupling constants and the field operators, 
{y\71q,gI"^} {x\71^,hI"^}. The full interaction Lagrangian is then given by adding the 
of both channels. 



6. Matching relations 

In the pion-nucleon sector, the coupling constants of the nonrelativistic Lagrangian, C, and 
dI*^^ can be expressed in terms of pion-nucleon scattering lengths of the S-wave and P-wave, «()+ 
and ai± and effective range parameters fto+> respectively. Adopting the notation of Ref. [1 1], in the 
isospin limit, the isospin decomposition of the nN scattering amplitudes reads 

Tnn+^nn+ ~ Tpn-^pn- — -\-T . (6.1) 

Defining .yV = 4n{mp+Mji), one finds 

Co = 2^aQ^, Ci=2V2^aQ^, C2 = 2 ^ (a^^ + a^^) , 

C3=Q), C4 = Ci, C5=C2. (6.2) 

(k) 

The matching conditions for the D] are given in a generic form only. The isospin index of the 
threshold parameters can be inferred from Eq. (^?I|)^ 

\2mpM„ 

Dp^ =2^ (ai_ -£?!+). (6.3) 

Here, higher order terms in the threshold parameters have been dropped. The corrections to these 
relations which appear due to isospin breaking have to be calculated within the underlying rela- 
tivistic theory. For the C,, they can be found in Refs. [12-15]. Note that the second line in Eq. ( ^!2| ) 
is only true in the isospin limit. 

(n) (n) 

The constants Gj and H- on the other hand are related to the threshold parameters of the 
electric and magnetic multipoles of the pertinent channel. In the isospin limit, the expansion of the 
real part of the multipole Xi± close to threshold is written in the form 

00 

ReX^±(^) = £X/±,2yt|q|'+^^ (6.4) 

k=0 



which defines the threshold parameters Xi±^2k- In the following, the relations of the coupling con- 
stants g|"' to these threshold parameters is given at leading order in the pion nucleon threshold 
parameters. Since the nonrelativistic theory is not suited for the study of the dependence of the 



Note that we use the Condon-Shortley phase convention. 
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multipoles on |k|, in this analysis, all vectors k are turned into unit vectors by the pertinent redefi- 
nition of the coupling constants, 

Cf^ = .Aok^"Gi , ^ = 4n{mp+M^o) . (6.5) 



Note that the higher order corrections due to Eq. ( |4.1[ ) have to be taken care of in the matching 
relations. Again, these relations pick up isospin breaking corrections which have to be evaluated in 
the underlying relativistic theory. 

Only the matching equations for the couplings of the Lagrangians if^'"'and^^*^^ are indicated 
here. The remaining relations can be found in Ref. [5]. To ease notation, Xi± = X,± o is used. 

Go = 2£o+, Gi =6(£+i+M+i), 

G2 = -2(M_i+2M+i), G3=6(£i+-Mi+). (6.6) 

For the coupling constants //, the algebraic form of the relations is identical. However, the multi- 



poles of the pertinent channels appear and the masses in Eq. (5.5) have to be adjusted. 



All coupling constants are assumed to be real. See Ref. [5] for a discussion of this issue. 
7. Results 

In the following, we provide the expressions for the electric and magnetic multipoles for 
Z = 0, 1 and Mi± for Z = 1 . The result is written in the form 

X/,±(.) =Xlr{s)+xl^''°^{s)+xf^°''^{s)--- (7.1) 

where s = {pi +k)'^ and the ellipsis denote higher order terms in the expansion in e and a. 

7.1 Tree-level 

The tree level result can be written in the form Xl^_^'^{s) = Xi_^_q' +X'i_^_ jQ^^' + • • • with the 
coefficients 

■^0+ = ^0 , 3£'o_|_ 2 = G4 — 3G5 + Ge — Gg , 

6Mi_|_ = Gi - G3 , ^{+,2 = - + i^Gio + 13G12 + gGi3 - ^Gi4 , 

3M[_ = G3 — Gi — 3G2 , ^[-,2 = 5^9 — 5G10 + G11 — |Gi2 



Gi3 + j^Gi4 , 



— Gi +G3 , ^{+,2 — + + j^Gn — \Gn 

+ 3^Gi4— ^Gi5. (7.2) 
One observes that D-waves appear naturally at order in this framework (see also Ref. [16]). 
7.2 One-loop 

All the one-loop contributions are proportional to the basic integral 



i{2n)D 2(0a {\)2(0b (P - 1) {coa (1) -h){cOh{V-\)-Po + k)' 
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Figure 1: One- and two loop topologies needed to calculate the amplitude. The double line generically 
denotes a nucleon, the dashed line a pion and the wiggly line indicates the external electromagnetic field. 



W±(p) = jM2+p2 



I = n,p 



In the limit D ^ 4, 

Jab{P'-) 



p2 = p2 _ p2 ^ 



^[s-{ma+M^tY){s-{ma-Mj,>,f), 



(7.3) 



(7.4) 



which is a quantity of order e. The one-loop result up to and including order O(ae^) reads 



(Pu Pn\ 
Pn Pii 
Psi P32 

\Ai P42 J 



Jn+is) 



(7.5) 



The elements P,^. are functions of the pion momentum q and the coupling constants of the La- 
grangian, 



Pu = GoCo + [CoE'^lJ - 2D}; 'Go) , 
Pu = CiHq + h\s, m,,M^,) (Ci^+ti'' - Df )//o ) - q'D\'>Ho , 



.2n{2), 



I8P22 = h^{s,m„M^,) (D['^-DfA (Hi-H^ 



18P4i=q' (G1+G3), 

I8P42 = h\s,m„M^,) (d['^-D^^A {Hi +H^) , 



(7.6) 



where E^^J^i denotes the pertinent coefficient of the tree level result of channel (c), see Eq. ([7^), 



As 



(7.7) 



and h^{s,mc,Mj^d) is given by 

h^{s,mc,Mj^d] 

which is a quantity of order e^. Eq. (7.5) and ( [7.6[ ) clearly show the advantage of the nonrelativistic 
description: The strength of the cusp at leading order is parameterized in terms of the coupling 
constant Ci and the ratio Hq/Gq. 
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8. Phase of it^p n^p scattering 

In an isospin symmetric world, the phase of the multipole £"0+ is directly related to the phase 
shift of the 5- wave of pion-nucleon scattering by virtue of the Fermi -Watson theorem [17]. In 
Ref. [7], it is shown with a coupled channel 5-matrix approach that to leading order in e, below the 
n~^n threshold, the phase of the S wave of Tr^jo n^p scattering is equal to the phase of £'o+> 



tan 5 0^ = tan ^^^^ = tan Se^^ ■ (8.1) 
Ke£,o+ 



The framework developed here allows one to test this statement order by order in the perturba- 
tive expansion. To this end, the phase of £0+ below the second threshold is calculated up to and 
including 0{a^e'^), 

tan 8Ea+ = Co Im/po + C?7„+Im7po - ID^^hmJpQC^ (8.2) 
- 2CiDf ^7„+Im7poq^ - 2CiDfhn+\mJpQh^{s,mn,M7t) + ■■■ ■ 



Calculating n^p —>■ n^p scattering to the same order with the Lagrangian given in Eq. (5.1), one 



finds that the phase of the 5- wave below the second threshold is indeed equal to Eq. (8.2). However, 
the main object of interest is the phase of the S'-wave of Ti^p —>■ n^p scattering in the isospin 
symmetry limit, 

tan Spjfi^p^ = ColmJpo + -^lmJ„+ ~ a^^q + 2-2±q + • • • , (8.3) 

Co ^0+ 

which does not agree with Spj^a^pj^o in the presence of isospin violations already at leading order. 



9. Summary and conclusion 

We study the photoproduction reaction of pions on the nucleon using a nonrelativistic frame- 
work. The electric and magnetic multipoles £"/+ for / = 0, 1 and Mi± are calculated in a systematic 
double expansion in the final state pion- and nucleon momenta (counted as a small quantity of 
order e) and the threshold parameters of nN scattering (denoted by a). Explicit representations for 
the multipoles up to and including and e^a are provided. The corresponding two-loop results as 
well as a expressions for the multipole amplitudes in the remaining three reaction channels can be 
found in Ref. [5]. 

The representation is valid in the low energy region, at least up to a photon energy in the lab 
frame of Ey = 165MeV. It accurately describes the cusp structure and allows one to determine the 
pion-nucleon threshold parameters from experimental data. 

The relation of the phase of the electric multipole £0+ in the (pO) channel to the phase of the 
S-wave of n^p —>■ n^p scattering is discussed in the presence of isospin violation. A relation found 
in earlier work [7] is confirmed. We stress that the relation does not allow one to obtain the phase 
of n^p —>■ 7r"/7 scattering in the isospin limit. 
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